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Abstract This paper presents a comprehensive analysis of the degree statis- 
tics in models for growing networks where new nodes enter one at a time 
and attach to one earlier node according to a stochastic rule. The models 
with uniform attachment, linear attachment (the Barabasi- Albert model), and 
generalized preferential attachment with initial attractiveness are successively 
considered. The main emphasis is on finite-size (i.e., finite-time) effects, which 
are shown to exhibit different behaviors in three regimes of the size-degree 
plane: stationary, finite-size scaling, large deviations. 
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1 Introduction 

Complex networks have attracted much attention over the last decades. They 
provide a natural setting to describe many phenomena in nature and society [TJ 
[21I31H1E] ■ One of the salient features of most networks, either natural and arti- 
ficial, is their scalefreeness. This term refers to the broad degree distribution 
exhibited by these networks. The probability that a node has degree k (i.e., 
is connected to exactly k other nodes) is commonly observed to fall off as a 
power law: 

fk~k-\ (1.1) 

This power-law behavior, which holds in the limit of an infinitely large network, 
will be referred to hereafter as 'stationary'. The exponent usually obeys 7 > 2, 
so that the mean degree of the infinite network is finite. Growing networks with 
a preferential attachment rule, such as the well-known Barabasi- Albert (BA) 
model [M7], have received a considerable interest, as they provide a natural 
explanation for the observed scalefreeness. The observation that preferential 
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attachment generates a power-law degree distribution actually dates back to 
much earlier works 019]. 

Scalefree networks, being chiefly characterized by the exponent 7 of their 
degree distribution, are therefore somewhat similar to equilibrium systems at 
their critical point. As a consequence, finite-size (i.e., finite-time) effects can 
be expected to yield important corrections to the asymptotic or stationary 
form (jl.ip of the degree distribution. These effects are one of the possible 
causes of the cutoff phenomenon which is often observed in the degree dis- 
tribution of real networks |10j . More precisely, the largest degree of a 
scalefree network at time n can be estimated by means of the following argu- 
ment of extreme value statistics: it is such that the stationary probability of 
having k > fc*(ri) is of order l/n. The largest degree thus grows as a power 
law [TOEl]: 

fc*(n)~n", f=— !-r. (1.2) 
7- 1 

This growth law is always subextensive, because one has 7 > 2, so that v < 1. 
The cases 2 < 7 < 3 (i.e., 1/2 < v < 1) and 7 > 3 (i.e., < v < 1/2) 
however correspond to qualitative differences, especially in the topology and 
in the various dimensions of the networks [12] . 

The goal of this article is to provide a systematic analysis of the degree 
statistics of growing network models at a large but finite time n. Both the 
age-resolved distribution fk(n, i) of the degree of node i at a later time n and 
the distribution fk{n) of an unspecified node at time n will be considered 
throughout. Several works have already been devoted to this problem, both 
for growing networks with preferential attachment [TT1IT31I141I151I161I17U18] and 
for related models of random graphs and other structures jT9][20] . The present 
work aims at being systematic in the following three respects: 
• Models. This work is focussed onto growing network models where a new node 
enters at each time step, so that nodes can be labeled by their birth date n, 
i.e., the time they enter the network. Node n attaches to a single earlier node 
(i = 1, . . . , n — 1) with probability p n ^. The attachment probabilities and the 
initial configuration entirely define the model. The network thus obtained has 
the topology of a tree. The degrees ki(n) of the nodes at time n obey the sum 
rule 

n 

^fci(n) =2L(n), (1.3) 

where L(n) is the number of links of the network at time n. 
We will successively consider the following models: 

- Uniform attachment (UA) (Section 2). The attachment probability is inde- 
pendent of the node, i.e., uniform over the network. This model is not scalefree. 
Its analysis serves as a warming up for that of the subsequent models. 

- Barabdsi- Albert (BA) model (Section 3). The attachment probability is pro- 
portional to the degree ki(n) of the earlier node. This well-known model [6l[7] 
is scalefree, with exponents 7 = 3 and v — 1/2. 
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Table 1.1 Various characteristics of the network for both initial conditions. The listed 
results hold irrespective of the attachment rule. 



Initial condition 



Case A 



Case B 



Topology 

Number of links at time n 
Mean degree at time n 



tree 

L( A )(n) =n-l 
<fc( A )(n)} = 2-2/n 



rooted tree 
L(B)(„) = n - 1/2 
(fc( B )(n)) =2-l/n 



Degrees at time 1 

and generating polynomials 







*i B) (l) = 



1 



Degrees at time 2 

and generating polynomials 



n(A) 



fc^ B) (2) = 2, 4 B) ( 2 ) = 1 
F 2 (B) (a;) = \x(x + l) 



- General preferential attachment (GPA) (Section 4). The attachment prob- 
ability is proportional to the sum fcj(n) + c of the degree of the earlier node 
and of an additive constant c > — 1. This parameter, representing the initial 
attractiveness of a node [11] . is relevant as it yields the continuously varying 
exponents 7 = c + 3 and v = l/(c+2). The BA and UA model are respectively 
recovered when c = and c — » 00. 

• Regimes. For each model, the following three regimes will be considered: 

- Stationary regime (k <C fc*(n)). The degree distribution is essentially given 
by its stationary form (|1.1[) . to be henceforth denoted by fk,st&t, in order to 
emphasize its belonging to the stationary regime. 

- Finite-size scaling regime (k ~ fc*(n)). In the scalefree cases, the degree 
distribution obeys a multiplicative finite-size scaling law of the form 

fk(n) w /fc.stat^f-^r) • (1.4) 

- Large- deviation regime (k+in) -C ~ n). The degree distribution is usually 
exponentially small in n. 

• Initial conditions. We will consider the following two initial conditions: 

- Case A. The first node appears at time n = 1 with degree fci(l) = 0. This 
prescription is natural because the first node initially has no connection. All 
subsequent nodes appear with degree k n (n) = 1. In particular, at time n = 2 
the second node connects to the first one, so that fci(2) = ^(2) = 1. The 
configuration thus obtained is the dimer configuration used e.g. in [TMll)] . At 
time n, the network has L(n) = n — 1 links. It has the topology of a tree. 

- Case B. The first node now appears at time n = 1 with degree fci(l) = 1. 
This formally amounts to saying that this node is connected to a root, which 
does not belong to the network. It is natural to associate half a link to this 
fictitious connection. At time n = 2 the second node connects to the first one, 
so that fci(2) = 2 and ^(2) = 1. At time n, the network has L(n) = n— 1/2 
links. It has the topology of a rooted tree. 

Table fTTTI summarizes various characteristics of the network for both initial 
conditions, whereas Figure [TTT1 illustrates the first three steps of the network 
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Fig. 1.1 First three steps of the construction of the network (upper panel) and corre- 
sponding interacting particle representation (lower panel) for both initial conditions. 



construction. The upper panel shows the networks with their nodes and links. 
The lower panel shows the corresponding representation as an interacting par- 
ticle system, where each node is viewed as a site occupied by a number of 
particles equal to its degree. The total number of particles in the system is 
therefore 2L(n). The information about the topology of the network, and es- 
pecially about the genealogy of the nodes, is lost in the interacting particle 
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representation, but this information will not be used in the present study which 
is focussed on the statistics of degrees. 



2 The uniform attachment (UA) model 

The uniform attachment (UA) model is the simplest of all: the attachment 
probability is chosen to be uniform over all existing nodes. This section is 
devoted to an analytical study of the distribution of the degree of a fixed node 
and of an unspecified node, exactly taking into account fluctuations, finite-time 
effects, and the influence of the initial condition. 



2.1 Degree statistics of a fixed node 

We start with the study of the distribution of the degree fcj(n) of node i at 
time n. The node appearing at time n > 2 links to any of the n — I earlier 
nodes (i = 1, . . . , n — 1) with uniform probability 

Pn,i = — ^r. (2.1) 
71 — 1 

If we define the degree increment of node j at a later time j > i as 

m = hu) - ki (j - 1) = { I ^ * probability *><> (2.2) 

the degree ki(n) of node i at a later time n is given by 

n 

h(n) = ki(i)+ ]T h(j), (2.3) 
j=i+i 

with fej(i) = 1, except for i = 1 in Case A, where fci(l) = (see Table [TTT|) . 
The mean degree (fcj(n)) therefore reads (i > 2) 

n -. 

(h(n)) = 1 + V ■— ■ = ff n _ a - Hi_ t + 1 w In - + 1, (2.4) 

7—1 7 

where the harmonic numbers H n are defined in (|2.20p . 

The distribution fk(n,i) — Prob{fcj(n) = k} can be encoded in the gener- 
ating polynomial 

n 

F n ,i(x) = (x k ^) =Y,fk{n,i)x k . (2.5) 
fe=i 

As a consequence of (|2.3I) . we have 



F nji (or)=^« [] (2.6) 
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where the characteristic function of the degree increment h{j) assumes the 
simple form 

(^W) = l + (x-l) Wii = £±iz2, (2.7) 
irrespective of i. We thus get (i > 2) 

F n M - f -Trfl"- 1 ? , (2-8) 

[n — 1)11 (x + % — 1) 
whereas only F n> i(x) depends on the initial condition according to 
( A) , , _ r(x + n-l) (B ) _ xr{x + n-l) 

[x) ~ (n-iy.rix)' n > l[ ) -(n-iy.r(x)- 1 yj 

Throughout the following, the superscripts (A) and (B) mark a result which 
holds for a prescribed initial condition (Case A or Case B). 

The product form (|2.6p implies that the generating polynomials of node i 
at times n and n + 1 obey the recursion 

F n+ i,i(x) = (a^" +1 >) ^(i) = a: + ^~ 1 F nii (a:). (2.10) 
The probabilities fk(n,i) therefore obey the recursion 

fk{n + l,i) = -f k -i{n,i)+(l--)f k (n,i), (2.11) 
n \ n J 

with initial conditions given in Table [TTT| i.e., 

fk(i,i)=S kt i (t>2), /^ A) (l ! l) = 4,o, /f ) (l,l) = 4,i. (2.12) 

The master equations (I2.11[) can be directly written down by means of a simple 
reasoning. They provide an alternative way of describing the evolution of the 
degree distribution of individual nodes. 

The degree distribution encoded in (|2.8|) has the following characteristics. 
The degree of node i at time n ranges from the minimal value 1 to the maximal 
value n + 1 — i. These extremal values occur with probabilities 

h{n,i) = — \, W,(v) = |^|. (2.13) 
n — 1 (n — 1)\ 

The mean and the variance of the degree can be obtained by expanding the 
result (|2.8p around x = 1, using 

(x K ) = l + (x-l)(K) + ^(x-l) 2 (K 2 -K) +•••, (2.14) 

var K+(K) 2 -{K) 

where K is any random variable taking positive integer values. We thus get 
(h(n)) = ff n _i-Fi_i + l, 

(2.15) 

varfci(n) - ff n _i - - H^i + h\% 
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where the harmonic numbers H n and Hn are defined in (|2.20p . The above re- 
sults hold irrespective of the initial condition. The first one coincides with (|2.4p . 

In the scaling regime where both times i and n are large and comparable, 
introducing the time ratio 

n 

z = ->\, (2.16) 

i 

the expressions (|2.15[) yield 

(ki(n)) s» hiz+ 1, varfcj(n) w In 2. ( 2 -17) 

In deriving the above results, we have used the asymptotic behavior of the 
digamma function W(x) — r'(x)/r(x) and of the trigamma function \P'(x) as 
x — > oo: 

W(x) = lnx - — H , ^'(sc) 

2a; 

as well as their values at integers: 

!P(n) = fln-i - 7*, *'{n) 

where 

= ^ 2) =E4 (2-20) 

i=i i=i 

are the harmonic numbers of the first and second kind, and 7 E is Euler's 
constant. 

The entire degree distribution can be characterized in the scaling regime. 
Equation (|2.8|) indeed yields 

F n>i (x) ^xe^" 1 ) 1 " 2 , (2.21) 

irrespective of the initial condition. We recognize the generating function of a 
Poissonian distribution with parameter A = hiz, up to a shift by one unit. We 
thus obtain [151I2T] 

Clnz^" 1 



= - + 5-j + ---, (2.18) 
= 4-^ (2-19) 



2.2 Degree statistics of the whole network 

We now turn to the degree distribution of the whole network at time n, fk(n) = 
Prob{fc(n) = fc}, where k(n) stands for the degree of an unspecified node. We 
have 

1 ™ 

h{n) = -Y j f k {n,i). (2.23) 

i=l 

The corresponding generating polynomials, 

n i n 

F n (x) = (x k( -^)=Y^fk{n)x k = -Y / F n 4x), (2.24) 

k=l i=l 



s 



obey the recursion 

(n + l)F n+1 (x) = (x + n-l)F n (x) + x, (2.25) 

or equivalently 

(n + l)/ fe (n + 1) = f k -i(n) + (n - l)f k (n) + 5 kA , (2.26) 
with initial conditions given in Table [TTTj i.e., 

jf>(l)=4,o, /i B) (l) = 4,i- (2-27) 
The recursion (|2.25[) has a non-polynomial solution, independent of n, 

F stat (x) = (2.28) 
2 — x 

describing the stationary degree distribution on an infinitely large network: 

/fc,stat = ^ (*>1). (2.29) 

The solution of (|2~25|) reads 

P ( A V^ - x 4. 2 ( l ~ x ) r(x + n-i) 

" ( ' - 2 - x + 2-x n\r(x) ' 
p (B), , _ x x(l - a) r(x + n - 1) 
" [) ~2-x + 2-x n\r(x) ■ 



(2.30) 



The polynomials Fn (x) and F n (x) have respective degrees n—1 and n. 
The first of them which are not listed in Table fTTT] read 



F 3 (A) (x) = \x{x + 2), FP(x) = \x(x 2 + 2x + 3), 

(2.31) 

F 4 (A) (x) = ^ x(x 2 + 4a; + 7), F 4 (B) (x) = £ x(x + 3)(x 2 + x + 4). 

The degree /c(n) at time n ranges from the minimal value 1 to the maximal 
value n—1 (Case A) or n (Case B). These extremal values occur with the 
following probabilities (n > 2) 

^w-5 + 55^ij- >*»4 J^w-a- ^»>-a- 

(2.32) 

We now turn to the finite-size scaling behavior of the degree distribution 
when both k and n are large. As anticipated in the Introduction, it is to be 
expected that the probabilities fk(n) are close to their limits {fk(n) w /fe, s tat) 
for ?i large at fixed degree k, and more generally in the stationary regime 
where k is much smaller than some characteristic crossover degree fc* (n) . Con- 
versely, the probabilities fk(n) are expected to be negligible (fk(n) <§C /fc.stat) 
for fc large enough at fixed time n, and more generally in the large-deviation 
regime where fc*(n) <C fc ~ n. The crossover scale fc*(n) can be estimated as 
fe*(n) » (ki(n)} (see (|2.4p ). Nodes with highest degrees are indeed typically 
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k/ln n 



Fig. 2.1 Plot of the ratios i?fc(n) against k/lnn (see 1 12, 340 . for the UA model with initial 
condition A, at times n = 10 3 (empty symbols) and n = 10 6 (full symbols). 



expected to be the oldest ones. An alternative route consists in using the ar- 
gument of extreme value statistics alluded to in the Introduction: the largest 
degree fc* at time n is such that the stationary probability of having k > k+ is 
of order 1/n. Both approaches consistently yield 

fc*(n)~lnn. (2.33) 
Finite-size effects are best revealed by considering the ratios 

Rk(n) = |^ - 2 k f k (n). (2.34) 

These ratios are expected to fall off to zero for k of the order of fc*(n) ~ Inn. 
Figure [27X1 shows a plot of the ratios Rk(n) against k/ Inn, for times n = 10 3 
and n = 10 6 in Case A. Numerically exact values of the /fc(n) are obtained by 
iterating (|2.26[) . A steeper and steeper crossover is clearly observed. 

In order to get some quantitative information on the observed crossover, 
it is advantageous to introduce the differences dk(n) — Rk-i(n) — Rk(n) for 
k > 2, completed by rfi(n) = 1 — i?i(n), i.e., Ro{n) — 1. Although the c4(n) 
are not positive, most of them are, and they sum up to unity, so that it is 
tempting to think of them as a narrow probability distribution living in the 
crossover region. The generating function of the dk(n) reads 

D n (x) = d k (n)x k = (x- l)F n (2x) + x. (2.35) 
fe>i 
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The above picture suggests to define the crossover scale as the first moment 

K = M(n) = kd k( n ) = Ki 1 ), (2-36) 

fc>i 

and the squared width of the crossover front as the variance 

a 2 (n) = k 2 d k (n) - M (n) 2 = D' 7 [{1) + /i(n) - M (n) 2 . (2.37) 

k>l 

Equations (l2~30]) . (j2~35]) yield 

p (A) (n) = 2H n m 2(lnn + 7E ), ^(n) = 2ff„ + U 2(ln fl + 7E ) + l, (2.38) 
and 

a 2 (n) = 2H n - 4H® » 2(lnra + 7e - tt 2 /3), (2.39) 

the latter result being independent of the initial condition. 

The crossover scale is thus k+ w 2 Inn, whereas the width of the crossover 
front grows as cr(n) ~ (2 Inn) 1 / 2 . These predictions are in agreement with the 
observations which can be made on Figure [2~Tl namely that the crossover takes 
place around k/hxn = 2, and that it becomes steeper at larger times, as its 
relative width falls off, albeit very slowly, as (mn)^ 1 / 2 . 

Another illustration of finite-size effects is provided by the complex zeros 
of the polynomials F n (x) . The location of these zeros indeed shows how fast 
the degree distribution of finite networks, encoded in the polynomials F n (x), 
converges to the stationary distribution, encoded in the function F sta t (%) ■ For 
n > 2, Fn (x) and F„ B \x) have one trivial zero at x — 0, and respectively 
n — 2 and n — 1 non-trivial ones. The explicit expressions (|2.30|) allow one 
to find the asymptotic locus of the zeros as follows. The most rapidly varying 
part of these results is the rightmost ratio, so that the zeros are asymptotically 
located on the curve with equation \r(x + n— l)/(n\r(x))\ = 1. Setting 

x = <, (2.40) 

and using Stirling's formula, we can recast the above estimate as 

Re[(l + £)M1+0-£1*£]=0. (2.41) 

The non-trivial zeros of the polynomials F n (x) are thus predicted to escape 
to infinity linearly with time n. Once rescaled by n according to (|2.40p . they 
accumulate onto a well-defined limiting curve in the complex £-plane. This 
curve, with equation (|2.41|) . has the shape of a lens connecting the points — 1 
and 0. Figure 12.21 illustrates this result with data at time n = 50 for both 
initial conditions. The polynomials F n (x) converge to the stationary function 
-Fstat(^) whenever the complex ratio £ = x/n lies within the lens. Otherwise 
they diverge exponentially with n. 
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Fig. 2.2 Plot of the non-trivial zeros of the polynomials F n (x) for the UA model, in the 
complex plane of the rescaled variable £ = x/n. Symbols: zeros for n = 50 in Case A (empty 
symbols) and Case B (full symbols). Line: limiting curve with equation 112.411 1. 



A related issue concerns the behavior of the probability (n) of having a 
very large degree, of order k ~ n, much larger than ~ Inn. Consider- 

ing Case A for definiteness, the expression (|2.30p leads to the exact contour- 
integral representation 

f (A) (n)-l dx (— I 2(1 ~ X) r (s + "-l) \ f2 42) 

The presence of gamma functions suggests to look for a saddle point x s pro- 
portional to n. Setting £ = k/n, we indeed find x s — n/v, where £ and v are 
related through 

hx(v + 1) , 
C = — -. 2.43 

v 

We thus obtain the following large-deviation estimate 

h (n) ~ exp (-n(C In n + S(())) , (2-44) 

where the exponent has a usual contribution in n and a less usual one in 
nlnn. The term linear in n involves a large-deviation function S(C), which is 
obtained in the following form, parametrized by v: 

5(C) = -(vlnv- \nv\n(v + l)-(v+ 1) lnfu + 1)). (2.45) 

v 

This function decreases from 5(0) — to 5(1) = — 1. The resulting behavior 
at £ = 1, i.e., exp(— n(lnn — 1)), is in agreement with the inverse factorial 
expressions (|2.32p . 
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3 Linear preferential attachment: the Barabasi-Albert (BA) model 

The Barabasi-Albert (BA) model is the simplest of the models with preferen- 
tial attachment: each new node connects to earlier nodes with a probability 
proportional to their degrees. The probability that node n connects to an 
earlier node i thus reads 

h{n- 1) 

Pn >< = zfr=T)> (3 - 1} 

where ki(n — 1) is the degree of node i at time n — 1, i.e., before node n enters 
the network. The partition function in the denominator, 



Z{n) = ^2 k i{n) = 2L(n) (3.2) 

i=l 

(see (|1.3p ). ensures that the attachment probabilities add up to unity. 

In the following we analyze the BA model along the lines of the previous 
section, keeping consistent notations as much as possible. The dependence of 
the attachment probability p n ^ on the degree ki(n — 1) however makes the 
problem more difficult than the previous one of a uniform attachment. 



3.1 Degree statistics of a fixed node 

Let us again begin with the distribution fk(n,i) = Prob{fci(n) = k} of the 
degree of node i at time n. 

A first estimate of the degree fcj (n) is provided by the following recursion 
relation for the mean degree (ki(n)), which is a consequence of (|2.2j) : 



(ki(n)) = (h(n - 1)) + (p^) = ^1 + z{n _ 1) ) (h(n - 1)). (3.3) 

In the scaling regime where both i and n are large, using the expressions of 
the partition function given in Table HTT1 i.e., 

Z {A) (n) =2n-2, Z {B \n) = 2n - 1, (3.4) 

the above relation becomes the differential equation 

d{h{n)) ^ {h{n)) 

dn ~ 2n ' ^' 5j 

which yields 

l) ■ (3 ' 6) 
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The generating polynomials F n j(x) and i^ n _|_i^(x) which encode the dis- 
tribution of the degree of node i at successive times n and n + 1 obey the 
recursion formula: 



= ((l + tx-lM^W) 



i.e. 



F w+M (s)=J' tM ( a ;) + V> , ' n f (3-8) 
Z(nj da; 

where Z(n) is given by Q3.4p . The probabilities fk(n,i) themselves therefore 
obey the recursion 

/ fe (n+l,i) = |=^/ fe _ 1 (n,i) + ^l-^y)/ fe (n,i), (3.9) 

with initial conditions (|2.12[) . The initial condition for Case A should be taken 
at time n = 2, in order to avoid indeterminate expressions, as Z^{\) = 0. 

In order to solve the recursion (|3.8|) . we perform the rational change of 
variable from x to u such that 

X U d d fn,n\ 

u = , x = , x(x — 1)— = — it— — . (3.10) 

1 — x u + 1 dx du 

Introducing the notation F Jly i(u) — F n ,i(x), the recursion (|3.8|) reads 

F n+U (u) = F nti (u) - " ■ (3.11) 
Z(nj du 

It is then advantageous to introduce the Mellin transform M n ^(s) of F n ^(u), 
defined as 

/>oo 

M n;i (s)= Fn.iMu-'-Mu. (3.12) 
Jo 

The inverse transform reads 

- / ^A^^(s)m s , (3.13) 

where C is a vertical contour in the complex s-plane whose position will be 
defined in a while. The virtue of the Mellin transformation is that the recur- 
sion l|3.1ip simplifies to 

M n+ x,i{s) =(l- M nti (s), (3.14) 
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with initial condition Mj^s) = Xq(s) for i > 2, with 

X (s)= x(u)u- s - 1 du= x- s (l-x) s - 1 dx = - (3.15) 

Jo Jo sm ns 

for < Res < 1. Hereafter the contour C is assumed to be in that strip. We 
thus get (i > 2) 

rW f ,_r(n-±-i)r(i-i) 



M$>{8) = , ' ' Xo{s), 

I (i — | — 1J 1 (n — 1) 

fBv , r(n-§-i) r(i-i) _ (3 - 16) 



These product formulas in the Mellin variable s are reminiscent of 

The mean and the variance of the degree of node i at time n can be 
extracted from these results as follows. The identity (|2 . 14[) yields 



F n ,i(u) = 1 1 — -5 1 (3-17) 



as u — > +00. Furthermore the coefficients of 1/u and 1/u 2 are respectively the 
residues of M rly i(s) at s = — 1 and s = —2. We thus obtain 

N 1 K " r{i-\)r{n-i) x 1 y " r(i)r(n-\) v ' 



and 

.(A) - . _ n - 1 



var fc^» = 2 — y - (*<(«)) ~ (*<>)>> 
varfcf >(n) =2^^ - (*f V)) 2 - (fcf } (n)>- 



(3.19) 



In the scaling regime where both times i and n are large and comparable, 
introducing the time ratio z = n/i (see (|2.16l) ). the above results yield 

(Jfcj(n)) w z 1/2 , varfc,(n) « z 1 / 2 ^ 1 / 2 - 1), (3.20) 

irrespective of the initial condition. The mean degree is in agreement with the 
estimate (|3.6[) . The entire degree distribution can actually be derived in the 
scaling regime. Equation (|3.16[) indeed yields 

M„,,( S )«z- s / 2 -^. (3.21) 
sm7rs 

We thus obtain 
and finally 

f k (n^)^z-^(l-z-^) k -\ (3.23) 

The degree distribution is therefore found to be asymptotically geometric, 
irrespective of the initial condition [T51I21) . 
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3.2 Degree statistics of the whole network 

We now turn to the degree distribution fk(n) = Prob{fc(n) = fc}, where k(n) 
stands for the degree of an unspecified node. 

The generating polynomials F n [x) obey the recursion 

(n + l)F n+1 (x) = nF n (x) + n ^zA^El + x (3.24) 

Z(n) ax 

where Z(n) is again given by (j3.4[) . and with initial conditions given in Ta- 
ble 11.11 The probabilities fk (n) themselves obey the recursion 

(n + l)f k (n + 1) = JLzl n f k -i(n) +(l- -^—) nf k (n) + 6 k>1 . (3.25) 



Z[n) J ^ ly ' \ Z{n) t 

The first generating polynomials which depend on the attachment rule read 

F ( 3 A \x) = \ x(x + 2), F^ B \x) = \x{2x 2 + 2x + 5), 

F 4 (A) (a;) = I x(x 2 + 2x + 5), F 4 (B) (x) = ^ x (6x 3 + 8x 2 + llx + 35). 

(3.26) 

The stationary degree distribution fk,st&t can be determined as the solution 
of (j3.25p which becomes independent of n for large n. We thus get 

(k + 2)A, stat = (k - 1)A- Mtat + 24,i, (3-27) 

hence [TTJ[T1[15] 

4 

/fe ' stat = k(k + i)(k + 2y (3 ' 28) 

An alternative approach consists in looking for the asymptotic generating func- 
tion F st nt(x) as the solution of (13.24[) which becomes independent of n for 
large n. We thus obtain the differential equation 

x(l - x)F^ t (x) + 2F stat (x) = 2x, (3.29) 

which has for solution 

2 2fl - r) 2 

W*)=3 ^-^-ln(l-a;). (3.30) 

x x 

Expanding this result as a power series in x allows one to recover (j3.28l) . 

The recursion (|3 . 24[) for the generating polynomials F n (x) can be solved 
along the lines of the above solution of the recursion (|3.8|) . The Mellin trans- 
forms M n {s) of the functions F n (u) = F n (x) obey the recursion 

(n + l)M„+i(s) = (l - -^-^j nM n (s) + X (s), (3.31) 
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with initial condition M% ( s ) — Afj (s) = Xq(s). Equation (|3.31[) has a 
special solution 

Z(n)X a (s) 

M n{s) = , : , — , 3.32) 
(s + 2)n 

whereas the general solution of the homogeneous equation shares the n-depen- 
dence of the expressions (|3.16l) . We thus obtain 




(3.33) 



(s + 2)n \ 

Mi*\s) = 1 ^L(2n-l + ( S + l i 
(s + 2jn V 

The common stationary limit of both expressions, 

M stat ( S ) = (3.34) 
s + 2 

is proportional to the special solution (|3.32[) . Recalling (|3.15[) . the inverse 
Mcllin transform of the above result, 

F rtat (u) = l-- + 4ln(« + l) 5 (3-35) 

is equivalent to (I3.30|) . 

The results p.33|) allow one to investigate, at least in principle, every fea- 
ture of the degree distribution fk (n) . Let us take the example of the probability 
/i(n) for a node to have degree one. The inverse formula (|3 . 13|) shows that 
this probability is equal to minus the residue of M n (s) at s = 1. The nature 
of the subleading corrections to the stationary value /i s tat = 2/3 depends on 
the initial condition. For Case A we obtain (n > 2) 

(A , _ 2(^-1) 4r(„ |) 2 2 + 4 + 
Sn 3y 7rrei (n — 1) 3 3n 3^/nn d '' 1 

More generally, all the probabilities /t(n) exhibit a singular correction in 
n -3/2_ (^ ase g ^he remarkable property that all the probabilities fk{n) 
are rational functions of time n. Their expansion at large n therefore only 
involves integer powers of 1/n. We have e.g. 

/x W ~ 3n "3 3n' 
(B ) = n 2 -2n + 3 = 1 _ J_ J_ (3 ' 37) 
2 l " j 3n(2n-3) 6 Yin 8n 2 

We now turn to the finite-size scaling behavior of the degree distribution 
when both k and n are large. The crossover scale k+(n) can again be esti- 
mated either using (|3.6|) or by the argument of extreme value statistics. Both 
approaches consistently yield 

K{n) ~ n 1 / 2 . (3.38) 
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We will now show that the degree distribution obeys the multiplicative finite- 
size scaling law 

/*(»)« /k.rtat <%), V=^J2, (3-39) 

where the scaling function ^(y) is non-universal, in the sense that it depends 
on the initial condition [TMTT] . The proof of the scaling behavior (|3.39[) and 
the determination of the scaling functions (y) and <&( B > (y) go as follows. 
Let us start with Case A. The second term of the expression Q3.33P for M„ A ^ (s) 
scales as a power law for large n: 

The inverse Mellin transform of the latter formula, 

i?(A) (u) w I f ^1 2( S + l)Jr„( S ) ( .1,2) s f3 41) 

describes the scaling behavior of Fn (u) in the regime where u and n are 
simultaneously large, with u/n 1 / 2 fixed. Finally, by inserting the above scaling 
estimate into the contour-integral representation 

permuting the order of integrals, opening up the it-contour and using 

du (u + i)*- 1 _ r(k) 



c 2tti u k - s+l r(s)r(k - s + 1) ' 

we obtain after some algebra the scaling form (|3 . 39[) . with 

1 - s\ /y\ s + 2 



(y ) = i + -= / — — -r 



(3.44) 



Case B can be dealt with along the same lines. We thus get the similar ex- 
pression 

^>- + / c ra4-i)(fr 

The above expressions can be evaluated by closing the contour to the right 
and summing the residues at the poles of the gamma functions. We thus get 



/tt ^ (2m + 3)m! V2 

<s(B)^ = 1 4. V (-l) m (m + l)(2m + 3) ( y_^ n+i 
yg> ^ (m + 2)\ 12/ 



(m + 2)! 

m>0 
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Fig. 3.1 Plot of the ratios f k (n)/f k . stat against the scaling variable y — fc/n 1 /^, for the BA 
model at time n = 10 3 (symbols) for both initial conditions. Full lines: asymptotic scaling 
functions ^^(y) and <p( B ' l (y). 



i.e., finally 



^ A )(y)=erfc(|) 



1 



y 



V/4 



V/4 



(3.47) 



where erfc denotes the complementary error function. The above expression 
for <2>( A ) can be found in [ITjUTT] , whereas that for <Z>( B ) has been shown in [T3] 
to hold for a slightly different attachment rule and initial condition. 

Figure |3~T1 shows a plot of the ratios fk(n)/fk,stat, against the scaling vari- 
able y = k/n 1 / 2 , at time n = 10 3 for both initial conditions. Exact values for 
the /fe(n) are obtained by iterating (|3.25|) . The data are well described by the 
predicted finite-size scaling functions $( A \y) and ^ B \y), shown as full lines. 

Both scaling functions share similar qualitative features. They start from 
the value 1 at y = 0, increase to a maximum, which is reached for y = 2 in 
Case A and for y = V6 in Case B, and fall off as exp(— ?/ 2 /4). They however 
differ at the quantitative level, both at small and large values of y: 



*<A)(y) = l + J* 

3V7T 



^(B)( y ) = 1+ V 

yy> 32 



2^ 



!L e -v 2 /\ 



(3.48) 



Apart from the additive constant 1, the scaling functions #( A ) and arc 
respectively an odd and an even function of y. This is the transcription in the 
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finite-size scaling regime of the phenomenon underlined when discussing (|3.36p 
and (|3.37|) . In particular, the first correction term at small y is in y 3 for <1>( A \ 
and in y 4 for <p( B \ 

Let us again close up with the location of the complex zeros of the poly- 
nomials F n (x). Considering Case A for definiteness, the result (|3.33p can be 
recast as the exact formula 

f5"M - — W») - i I ± JK*: 1 ^ (3-49) 



c 



i s + 2 r(l-f)r(n-l) sinks' 



The growth of this expression with n for a fixed value of the complex variable u 
can be investigated by means of the saddle-point approximation. The presence 
of gamma functions again suggest to look for a saddle point s s proportional 
to n. Skipping details, let us mention that we find s s « 2n/(l — it 2 ), so that 
the right-hand side of (|3.49p can be estimated as 

F n , sins (u) ~(x-^J , (3.50) 

with exponential accuracy. The asymptotic locus of the complex zeros is then 
naturally given by the condition that the above estimate neither falls off nor 
grows exponentially. We thus obtain |l — l/u 2 | = 1. The relevant part of this 
locus can be parametrized by an angle < 6 < 2ir as 

u = {1 - c- w y 1/2 , x= m . (3.51) 

This closed curve in the cc-plane has a cusp at the point x = 1, corresponding 
to the scaling regime, with a right opening angle. We have indeed x — 1 w 
(e 17r/2 6») 1/2 as 9 -> 0. Figure l3~2l illustrates this result with data at time n = 50 
for both initial conditions. The polynomials F n (x) converge to the stationary 
series -F s tat(^) whenever the complex variable x lies within the closed curve 
shown on the figure. Otherwise they diverge exponentially with n. 

The exponential estimate (|3 . 50[) has another virtue. By inserting it into 
the contour- integral representation (|3.42p , we obtain 

du (u+ l\ k f 1 

27T1 \ U J \ U- 



— ) (1-4) • (3.52) 



This integral can in turn be investigated by means of the saddle-point approx- 
imation. The result is the following large-deviation estimate 

f k {n) ~ exp(-nS(C)), (3.53) 
where C, = k/n, and where the large-deviation function reads 

5(C) - (1 - C) ln(l - C) - (2 - C) In (3.54) 
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The formula (|3.53[) describes, with exponential accuracy, the degree distribu- 
tion in the whole large-deviation regime where k and n are comparable. The 
quadratic growth S(() ~ C 2 /4 at small £ matches the fall-off of the finite-size 
scaling functions $ (A) (y) ~ $^(y) ~ exp(-y 2 /4) (see (|3.48jl ). The maximal 
value S(l) = In 2 describes the fall-off fk(n) ~ 2~" of the probability of having 
a degree k equal to its maximal value (k — n or k = n — 1). 



4 The general preferential attachment (GPA) model 

We now consider the general preferential attachment (GPA) rule, where the 
attachment probability to a node is proportional to the sum ki (n) + c of the 
degree of the earlier node and of an additive constant c, representing the initial 
attractiveness of the node [11] . This attachment rule interpolates between the 
uniform attachment rule, which is recovered in the c — > oo limit, and the BA 
model, which corresponds to c = 0. It can actually be continued on the other 
side of the BA model, as c can be chosen in the range — 1 < c < oo. The GPA 
model thus defined is scalefree for any finite value of c, with the continuously 
varying exponents 7 = c + 3 and v = l/(c+ 2). 

The probability that node n connects to an earlier node i thus reads 



_ ki(n — 1) + c 
'* ~ Z(n-1) 



(4.1) 
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where fc,(n — 1) is the degree of node i at time n — 1, and the partition function 
in the denominator, 

n 

Z(n) = ^2(h(n) +c) = 2L(n) + en (4.2) 

i=i 

(see (|1.3p ). ensures that the attachment probabilities add up to unity. 

In the following we analyze the GPA model along the very lines of the 
previous section. 



4.1 Degree statistics of a fixed node 

Let us again begin with the distribution fk{n,i) — Prob{fci(n) = fc} of the 
degree of node i at time n. 

A first estimate of the degree fa (n) is provided by the product formula (|3.3p 
for the mean degree (fci(n)), which still holds in the present case. In the scaling 
regime where both i and n are large, the latter relation becomes the differential 
equation 

d(h(n)) _ (kj(n))+c 

dn ~ (c + 2)n ' 1 ' 1 

which yields 

/T7 \ l/(c+2) 

(*i(n))«(c +1)(-) -c. (4.4) 

As anticipated, this expressions exhibits a power-law growth with exponent 
v = l/(c + 2) in the range < v < 1. 

The generating polynomials F n ^{x) and Fyj+i^a;) associated with the de- 
gree of node i at successive times n and n + 1 obey the recursion formula: 

F n+1>i (x) = = ^(»+D x Mn)\> 

= ((1 + (X-1K+I,i)z fci(n) ) 

= ((l + ^(fc^)+ C ))^ (n) ), (4-5) 



i.e. 



where 



Z^(n) = (c + 2)n-2, Z< B ' (n) = (c + 2)n - 1. (4.7) 
The probabilities fk{n,i) themselves therefore obey the recursion 

fk(n+l,i) = k+ z C { ~ ) 1 h-i(n,i) +(l- |^y) fk(n,i), (4.8) 
with initial conditions (|2.12[) . 
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In order to solve the recursion (|4.6jl . we again perform the change of vari- 
able (|3.10p from x to u, and set 

F n ,i{x) = (l-x)- c F n 4u). (4.9) 
The recursion (|4.6[) then reads 

F n+hi (u) = F n ,i(u) - -=^r [cfn, i (u) + u^^ J . (4.10) 

We then again introduce the Mellin transform M„ i j(s) of F ny i(u), so that the 
recursion (|4.10[) simplifies to 

M n+1>i (s) = (l - M n ,i{s), (4.11) 

with initial condition M,- ( j(s) = X c (s) for i > 2, with 

*.w = f *-d - *r +c - x = r(1 :; s) f;+ c) (4.12) 

Jo r ( c + 1 ) 

for — c < Re s < 1. Hereafter the contour C is assumed to be in that strip. We 
thus get (i > 2) 

m%(s)= ; + *J } c+ 2 2 { x c (s), 



r 



, I , (4-13) 

These product formulas are a generalization of (|3.16p . The mean and the 
variance of the degree of node i at time n can be extracted from these results 
as follows. The identity (|2.14[) now yields 

~ _ 1 (h(n)) + c (hjn) 2 ) + (2c + l)(h(n)) + c(c + 1) 

as u -» +oo. Furthermore the coefficients of this expansion are respectively 
the residues of M n: i(s) at s = — c, s = — c — 1 and s = — c — 2. We thus obtain 

— ) I {i 2 



(4.15) 



n(A), \i , -,n V c+2 y I c+2 

>)) = (c+1) 



r 



(fcf ) (n)) = ( C +l) 



r(n)r(»--i 5 
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and 



r(n)rli- -|g 
var^ (A) (n) = (c + l)(c + 2) ^ — 



-(fcrH> 2 -(2 C +i)(fcr(")>-c( C +i), 

varfc| B) (n) = (c+l)(c+2)- 



(4.16) 



r 



-<fcf ) (n)) 2 ~(2c + l)<A ;i (B) (n))-c(c+l). 

In the scaling regime where both times i and n are large and comparable, 
introducing the time ratio z — n/i (see (12.161) ). the above results yield 

(ki(n)) pa (c+l)z 1/(c+2) -c, varfci(n) « (c+1)z 1/(c+2) (z 1/(c+2) -1), (4.17) 

irrespective of the initial condition. The mean degree is in agreement with the 
estimate (|4.4[) . The entire degree distribution can actually be derived in the 
scaling regime. Equation (|4.13|) indeed yields 

M„,,( S ) « z-(*+ c )/( c + 2 ) X c (s). (4.18) 

We thus obtain after some algebra 



» -+T (4-19) 



and finally 



f k (n,i) « z-( c +D/( c + 2 )(l - z -i/(c+2)Nfc-i r ( fc + c ) (4.20) 

J v ' ; v 1 r(k)r{c + i) y ' 

This result allows one to recover both the Poissonian law (|2.22p in the c — > oo 
limit and the geometric one (|3.23|) as c = 0. 



4.2 Degree statistics of the whole network 

We now turn to the degree distribution of the whole network at time n, fk(n) = 
Prob{fc(n) = k}, where k(n) stands for the degree of an unspecified node. 
The generating polynomials F n (x) obey the recursion 

(n + l)F n+1 (x) = nF n (x) + "^"^ (cF n (x) + x^^) + x, (4.21) 



where Z(ri) is given by (|4.7p . and with initial conditions given in Table ITTTI 
The probabilities ,fk{n) themselves obey the recursion 

(n + l)/ fc (n + 1) = fc ^, C 7W fc -i(n) + f 1 - ^) n/ fc (n) + (4-22) 
Z(n) V z ( n )J 
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The first generating polynomials which depend on the attachment rule read 



Ft\x) = 


| x(x -t 


-2), 


F 3 (B) (x) = 


l 


x((c 


3(2c+3) 




1 


x ((c 


4(3c+4) 


F 4 (B) (x) = 


1 




4(2c+3)(3c+5) 



(4.23) 



+ (c + l)(7c + ll).x + (3c + 5)(4c + 7)). 

The stationary degree distribution fk, s tat can be determined as the solution 
of (|4.22[) which becomes independent of n for large n. We thus get 

(fc + 2c+2)/ fc , stat = (k + c - l)/ fc -i, st at + (c + 2)4,1, (4.24) 

hence [TTJIH] 

/ fcsta t= (C ! 2)r(2C + 3)r(fc + C) - (4.25) 
This result has a power-law decay at large k: 

A ,„ , <£±ffi±iU-<=«, (4 . 26) 
r(c + i) 

An alternative approach consists in looking for the generating function -F s t a t (x) 
as the stationary solution of (|4.2ip . We thus obtain the differential equation 

x(l - x)F^{x) + (2c + 2 - cx)F stat {x) = (c + 2)z, (4.27) 

which is equivalent to (14.241) . The solution 



2c+2 7o (i-y) c+3 

can be recast in terms of a hypergeometric function, which boils down to 
elementary functions whenever 2c is an integer. 

Throughout the regime where the degree k and the parameter c are both 
large and comparable, the expression (|4.25|) assumes a stationary large-devia- 
tion form, 

/fc.stat ~ exp(-c0(«;)), (4.29) 

where k = k/c, and with 

0(k) = (k + 2) ln( K + 2) - (k + 1) ln(/c + 1) - 2 In 2. (4.30) 

The linear behavior 4>(n) ~ Kin 2 as k — > matches the exponential de- 
cay (|2.29[) of the stationary distribution in the UA model, formally corre- 
sponding to c -> oo, whereas the logarithmic growth 4>{n) ps InK + 1 — 2 In 2 
as k — ► oo matches the power- law decay (|4.26[) . 
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The moments of the stationary distribution, 



k>l 



can be derived from (|4.24j) . which yields the recursion 

p-2 



(c + 2 — p)m p = c + 2 +pcm p _i + ^ ( P j (m 9+ i + cm,). (4.32) 

9=0 



We thus get 

2(3c + 2) 



mo = 1, mi — 2, m2 



c 



2(13c 2 + 17c + 6) 2(3c + 2)(25c 2 + 33c + 14) 
m3 = ^1) ' W4 = c(c-l)( C -2) ■ 



(4.33) 



and so on. The power-law decay (|4.26|) implies that the moment m p is conver- 
gent for c > p — 2. 

The recursion (|4.2ip for the generating polynomials F n (x) can again be 
exactly solved for a finite time n. The Mellin transforms M n (s) of the functions 
F n {u) = (l-x) c F n (x) obey 

(n + l)M n+1 (s) =(l- nM n (s) + X c (s), (4.34) 

with initial condition M^ A) (s) = m[ B \s) = X c (s). Equation (|43^1) has a 
special solution 

Z(re)X c (s) 
(s + 2c + 2)n 

whereas the general solution of the homogeneous equation shares the n-depen- 
dence of the expressions (|4.13l) . We thus get 



{s + 2c + 2)n 
x (c + 2)n-2 + 2(s + c+l) 



r ( i - —A fin 2 



c+2 / V c + 2 



w (s + 2c + 2)ti 

^ n _ s+c+l \ pi C I 

x I (c + 2)n - 1 + (s + c + 1) 



(4.36) 



26 



In order to illustrate these general results, let us again consider the proba- 
bility fi(n) for a node to have degree one. This probability is minus the residue 
of M n (s) at s = 1. For Case A we obtain (n > 2) 



/i (A) 0») = ro I (c + 2)n-2 + 2(c + 2) 




r(n-m) r {m) 

2(c + 2)r( 



r 

2c+2 
c+2 



r 



2e+3 

n c + 2 +•••), (4.37) 



'(a) 

whereas for Case B we obtain (n > 2) 

This rational expression for f\ (n) is however an exception. The probabilities 
fk(n>) indeed generically have a singular correction in n~( 2c+3 )/( c+2 ) for both 
initial conditions, whereas only /i (ft) and f^\n) are rational functions of 
time n. 

We now turn to the finite-size scaling behavior of the degree distribution 
when both k and n are large. The crossover scale fc*(n) can again be esti- 
mated either using (14. 4|) or by the argument of extreme value statistics. Both 
approaches consistently yield 

fc*(n)~n 1 /^. (4.39) 
The degree distribution obeys a finite-size scaling law of the form 

k 

fk(n) « A,stat sP(tf), y = 1/(c+2) , ( 4 - 4 0) 

where the scaling function <£(?/) again depends on the initial condition [T6||17j. 
The determination of the scaling functions (y) and <p( B ^ (y) closely follows 
the steps of Section 3.2. We thus obtain 



0(A) (y) = I 



f ds s + c+1 r(i- s ) +2c+2 

y j 



(c+ 2)/Y2c + 3) Jc 2?ri s + 2c + 2 r A, s_ 

1 c+2 



r 



ds s + c+1 r(l-s) 



yyj (c + 2)r(2c + 3) 7c 2tti s + 2c + 2 r fi^£ X 



c+2 



(4.41) 
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By closing the contours to the right, we can derive the following convergent 
series: 



2T 



2c+2 



$W( V ) = 1 + --\ c+2 ) v 2c +3 Y (rn + c+2)(-yr 

(c + 2)r(2c + 3) y 4'o(m + 2c+3)m!r(l-^ 

*(B) ( v ) = i + „ 2c+ 3 V (m + c + 2)(-»r 



(c + 2)r(2c+3) ^i( m + 2c + 3)m!r(-^)' 

(4.42) 

The above expression for <£( A ) can be found in [T7j , albeit not in a fully explicit 
form. It is also worth mentioning that the finite-size scaling function derived 
in [20j for asymmetric growing networks is different from the above one for 
generic values of the exponent v = l/(c + 2), although it coincides for v = 1/2 
with our result (|3T47|) for 

The expressions (|4.42[) suggest that the derivatives (y) and <p( B ^ (y) 
are somewhat simpler than the functions themselves. The factor (m + 2c + 3) is 
indeed chased away from the denominators under differentiation. The resulting 
series can be resummed by means of the identities 



(-y) m .. r ^ (4 ' 43) 



y sua = y 



c 27rl 



which are known e.g. in the theory of Levy stable laws. We are thus left with 
the following alternative contour-integral expressions for the derivatives: 



2F 



2c+2 



Both scaling functions start increasing from the value 1 according to the 
power laws 

2r( 2 ^i) 

<? (A) (y) = i + Vc+ 7 s y 2c+3 + ■■■, 

r( 2c + 4)r(f±i (4.45) 



r(2c + 3) J c 2m (4 44) 



c+1 



W + 2(c+2)3r(2c + 3) y ' 
go through a maximum, and fall off superexponentially as 

<pW(y)^2(c+2)Cr^)*(y), $M(y)nCr^)y#(y), (4.46) 
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with 



and 



= y 



2c+3- 



exp -(c + 1) 



c+2 ' 
c+1 



2c+3 



(2tt(c + l)) 1/2 (c + 2)2( c +i) r(2c + 3) 



(4.47) 



(4.48) 



Figure |4~T1 shows a plot of the scaling functions <p( A ^(y) and ^ B )(y) for (1) 
c = -1/2, i.e., t/ = 2/3; (2) c = 0, i.e., u = l/2 (the BA model); and (3) c = 1, 
i.e., i/ = 1/3. The figure demonstrates that the scaling functions present a high 
and narrow maximum for the smaller values of c, and a direct crossover from 
1 to for the larger values of c. These observations can be made quantitative 
by means of the pseudo-moments 



MP = 



POO poo 

/ <P'(y)y"dy^ P dy. 
Jo Jo 



(4.49) 



The integral formulas (|4.44|) allow one to evaluate these quantities explicitly: 



(A) 
Mp = 



(B) 
Mp 



(p + c + i)r(^M)r( P + 2 C + 3) 
(c + i)r(£±^)r(2 C + 3) 



(4.50) 



r(£±f±±)r(2c + 3) 
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• For large values of c (i.e., c — > oo), the model is close to the UA model. 
The analysis of the scaling functions will follow that of the ratios Rk (n) in the 
UA model, performed in Section 2.2. The crossover value of y, at which the 
functions exhibit a relatively sharp crossover from 1 to 0, can be estimated 
as Hi, i.e., 

^ = 2c + 2 7e + 2 + • ■ • , /4 B) = 2c + 2 7e + 3 + • • • , (4.51) 

which grows as 2c, irrespective of the initial condition. Similarly, the squared 
width of the crossover region can be estimated as the pseudo-variance a 1 = 
[i 2 - Mi, i-e-, 

ct 2(a) = 2c+4 7E + 2-27r 2 /3 + ---, cr 2(B) = 2c + 47 E + 3-2vr 2 /3 + - • • , (4.52) 

which also grows as 2c, irrespective of the initial condition. 

• For small values of c (i.e., c — > —1), the scaling functions exhibit a high 
and narrow peak around y = 1 . The position of the peak can be estimated as 
(y) = jti2/(2/ii), i.e., setting c = -1 + e, 

(y)^ = 1 + (3/2- 7e ) £ + ..., (j / )(b) = i + (2 - 7b ) £ + ..., (4.53) 

whereas the squared width of the peak can be estimated as vary = (4/Lii/i3 — 
3 M !)/(12 M ?), i.e., 

(A) 5 19 - 10 7e - 7T 2 2 
var y< A > = - e + e 2 + • • • , 

9 99 J ^ (4-54) 
rni * ZZ — o 7e — 7T 9 

var?/ B ) = - e + p e 2 + ■ ■ ■ , 

3 6 

and finally the area under the peak scales as fi±, i.e., 

/4 A) = - + 2 - 7e + • • • , /4 B) = - + 3 - 7e + • • • (4.55) 

e e 

We are thus left with the picture of a narrow peak around y = 1, whose width 
shrinks as e 1 / 2 and whose height grows as e -3 / 2 . 

Let us close up this section with the location of the complex zeros of the 
polynomials F n (x). The derivation of the estimate (|3.50p can be generalized 
to the present situation for arbitrary values of c. We are thus left with 

Fn«W~(i-r4^) , (4.56) 

again with exponential accuracy. The asymptotic locus of the complex zeros 
is therefore given by the condition |l — l/(— it) c+2 | = 1. The relevant part of 
this locus can be parametrized by an angle < 8 < 2ir as 

u=(l- e - w V mc+2) x = - (4 57) 

This closed curve in the cc-plane has a cusp at the point x = 1, corresponding 
to the scaling regime, with an opening angle equal to 7r/(c+2). We have indeed 
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-7 -5 -3-11 3 
Re x 

Fig. 4.2 Plot of the non-trivial zeros of the polynomials F n (x) in the complex x-plane. 
Symbols: zeros for n = 50 in Case A (empty symbols) and Case B (full symbols). Lines: 
limiting curves with equation 114.571 1. From the inside to the outside: c = (the BA model, 
already shown in Figure [3721 ) . c = 1 and c = 2. 



x-Iki (W 2 6>) 1/(c+2) as 9 -> 0. Figure l42l illustrates this result with data at 
time n = 50 for three values of c and both initial conditions. 

The exponential estimate (|4.56j) can again be recast into a large-deviation 
estimate for the probabilities fk(n) in the regime k ~ n, of the form 

A(n)~exp(-n£(C)), (4.58) 

with £ = k/n. The large-deviation function S(£) is obtained in parametric 
form: 

(c+2)(v-l) 

^ „c+2 _ 1 ' 

V 1 (4.59) 
S = ln(v c + 2 - 1) °~t 2 Uv - 1) ln(w - 1) + (v c+1 - l)v Inv) , 

where the parameter v in the range 1 < v < oo is the opposite of the saddle- 
point value of u in the contour integral generalizing (|3.52[) . 
The power-law behavior 

c+2 
( \ c+1 



5(C)«(c+l)l^l (4.60) 

at small C (corresponding to u — > oo) exactly matches the superexponential 
decay (|4.46[) . (|4.47[) of the finite-size scaling functions ^ A ^(y) and $( B \y). 
The maximal value 

S(l) =ln(c+2), (4.61) 
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corresponding to v — > 1, describes the exponential decay fk(n) ~ (c + 2) _ ™ 
of the probability of having a degree fc equal to its maximal value (k = n or 
fc = n — 1). 

5 Discussion 

In this paper we have presented a comprehensive study of finite-size (i.e., finite- 
time) effects on the degree statistics in growing networks. We have considered 
models defined by stochastic attachment rules, where nodes enter the network 
one at a time and attach to one single earlier node, so that the network has 
the topology of a tree. The present study thus generalizes and extends many 
results of References [TT]|I5P11[T51[TdT[T7|[T5] . 

We have successively investigated the uniform attachment rule (UA), the 
linear attachment rule of the Barabasi- Albert (BA) model, and a general 
preferential attachment rule (GPA) characterized by a continuous parame- 
ter c > — 1, representing the initial attractiveness of a node. The UA and BA 
models are recovered as two special cases, respectively corresponding to c — > oo 
and c = 0. The model is scalefree for any finite value of c, with the continuously 
varying exponents 7 = c + 3 and v = l/(c + 2). The continuous dependence of 
exponents on the parameter c, and the dependence of finite-size scaling func- 
tions on the initial condition (Case A or Case B in the present study), are 
two illustrations of the lack of universality which altogether characterizes the 
scaling behavior of growing networks. 

The GPA rule is actually the most general one for which the partition func- 
tion Z(n) (see (|4.2[) ) is deterministic, i.e., independent of the history of the 
network. Whenever the attachment probability has a non-linear dependence 
on the degree k, the partition function becomes a history-dependent fluctu- 
ating quantity, so that the analysis of size effects becomes far more difficult. 
The general case of an arbitrary attachment rule, growing either less or more 
rapidly than linearly with the degree, has been considered in several works [141 
[151H8] . Whenever the degree dependence of the attachment rule is asymptoti- 
cally linear, the resulting network is generically scalefree. The determination of 
the degree exponent 7 is however a highly non-trivial task in general (see [141 
15 for an explicit example). 

The present study has underlined the key role played by the typical value 
fe*(n) of the largest degree in a finite network at time n. In the UA model, 
ki,{n) grows logarithmically with time n. The situation is more interesting 
in the scalefree case, i.e., for c finite. The largest degree fc*(n.) grows as a 
subextensive power law with exponent v, and demarcates three regimes in 
the size-degree plane, where finite-size (i.e., finite-time) effects on the degree 
distribution fk(ji) have different forms. 

In the stationary regime {k -C fc*(n)), the degree distribution is very close 
to the stationary one, fk. ,stat - 

In the finite-size scaling regime (k ~ fc*(n)), the degree distribution obeys 
a multiplicative finite-size scaling law. As already noticed in several earlier 
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works, the finite-size scaling function <P depends on the initial condition im- 
posed on the network. This lack of universality holds for all finite values of 
the parameter c. Another feature of the finite-size scaling function is that it 
increases from its initial value <?(0) = 1, reaches a maximum, and stays above 
unity for a range of values of its argument y = k/k*{n), before it eventually 
falls off to zero. This non-monotonic overshooting behavior is however not 
mandatory. In this respect it is worth recalling the example of the so-called 
zeta urn model [2"2 ] l23 1 l24] . This mean-field interacting particle system with 
multiple occupancies possesses a continuous condensation transition at a fi- 
nite critical density. Its behavior right at the critical density shares a high 
amount of similarity with the present problem, including a power-law station- 
ary distribution with a continuously varying exponent, and finite-time scaling. 
The same results have been shown to apply to the dynamics of condensation 
in the zero-range process (ZRP) [25]. In the critical zeta urn and ZRP models, 
the finite-size scaling function is a monotonically decreasing function, so that 
<£(y) < 1 for all y. This does not contradict the conservation of probability: the 
excess probability is carried by smaller values of k, pertaining to the stationary 
regime. 

- In the large-deviation regime <C k ~ n), the degree distribution 

falls off exponentially in n. At variance with the finite-size scaling law, the 
corresponding large-deviation function is independent of the initial condition. 
The analysis of this regime has been shown to be closely related to the locus 
of the complex zeros of the generating polynomials F n (x), which have played 
a central role throughout this work. 

To close up, it is to be hoped that some of the concepts and methods used 
in the present work can be used to shed some new light either to other ob- 
servables in the network models considered here, such as e.g. the statistics of 
leaders and lead changes [21] , or to the degree statistics in more complex net- 
work models, such as e.g. the Bianconi-Barabasi (BB) model [26|i27j. where 
attachment rules involve the competing effects of dynamical variables (the 
node degrees) and quenched disordered ones (the node fitnesses). Depending 
on the a priori distribution of the random fitnesses, the BB model may possess 
a low-temperature condensed phase. Some features of the dynamics of the con- 
densed phase have been investigated recently, both at zero temperature [28], 
where the model is intimately related to the statistics of records, and at finite 
temperature [29] . 
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